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Soliton solutions of the KP equation have been studied since 1970, 
when Kadomtsev and Petviashvili proposed a two-dimensional non- 
linear dispersive wave equation now known as the KP equation. It 
is well-known that the Wronskian approach to the KP equation pro- 
vides a method to construct soliton solutions. The regular soliton 
solutions that one obtains in this way come from points of the totally 
non-negative part of the Grassmannian. In this paper we explain how 
the theory of total positivity and cluster algebras provides a frame- 
work for understanding these soliton solutions to the KP equation. 
We then use this framework to give an explicit construction of cer- 
tain soliton contour graphs, and solve the inverse problem for soliton 
solutions coming from the totally positive part of the Grassmannian. 

Introduction 

The KP equation, introduced in 1970 [1 , is considered to be 
a prototype of an integrable nonlinear dispersive wave equa- 
tion with two spatial dimensions. Concretely, solutions to 
this equation provide a close approximation to the behavior 
of shallow water waves, such as beach waves. Given a point A 
in the real Grassmannian, one can construct a solution to the 
KP equation [2]; this solution ua{x, y, t) is written in terms of 
a T- function, which is a sum of exponentials. More recently, 
several authors [3l [H [5l [6j have focused on understanding the 
regular soliton solutions that one obtains in this way: these 
come from points of the totally non-negative part of the Grass- 
mannian. 

The classical theory of total positivity concerns square 
matrices in which all minors are positive. This theory was 
pioneered in the 1930's by Gantmacher, Krein, and Schoen- 
berg, and subsequently generalized in the 1990's by Lusztig 
[II [81 , who in particular introduced the totally positive and 
non-negative parts of real partial flag varieties. 

One of the most important partial flag varieties is the 
Grassmannian. Postnikov [9 investigated the totally non- 
negative part of the Grassmannian (G'rfcn)>o, which can be 
defined as the subset of the real Grassmannian where all 
Pliicker coordinates are non-negative. Specifying which mi- 
nors are strictly positive and which are zero gives a decompo- 
sition into positroid cells. Postnikov introduced a variety of 
combinatorial objects, including decorated permutations, J- 
diagrams, plabic graphs, and Grassmann necklaces, in order 
to index the cells and describe their properties. 

In this paper we develop a tight connection between the 
theory of total positivity for the Grassmannian and the behav- 
ior of the corresponding soliton solutions to the KP equation. 
To understand a soliton solution ua{x, y, t), one fixes the time 
t, and plots the points where ua{x, y, t) has a local maximum. 
This gives rise to a tropical curve in the xy-plane; concretely, 
this shows the positions in the plane where the correspond- 
ing wave has a peak. The decorated permutation indexing the 
cell containing A determines the asymptotic behavior of the 
soliton solution ai y ^ ±oo. When t is sufficiently small, we 
can predict the combinatorial structure of this tropical curve 
using the d- diagram indexing the cell containing A. When 
A comes from a totally positive Schubert cell, we show that 
generically this tropical curve is a realization of one of Post- 
nikov's reduced plabic graphs. Furthermore, if we label each 
region of the complement of the tropical curve with the dom- 
inant exponential in the r-function, then the labels of the un- 



bounded regions form the Grassmann necklace indexing the 
cell containing A. Finally, when A belongs to the totally pos- 
itive Grassmannian, we show that the dominant exponentials 
labeling regions of the tropical curve form a cluster for the 
cluster algebra of the Grassmannian. Letting t vary, one may 
observe cluster transformations. 

These previously undescribed connections between KP 
solitons, cluster algebras, and total positivity promise to be 
very powerful. For example, using some machinery from total 
positivity and cluster algebras, we solve the inverse problem 
for soliton solutions from the totally positive Grassmannian. 

Total positivity for the Grassmannian 

The real Grassmannian Grkn is the space of all /c-dimensional 
subspaces of R^. An element of Grkn can be represented by 
a full-rank k xn matrix modulo left multiplication by nonsin- 
gular k X k matrices. 

Let (^^^) be the set of /c-element subsets of [n] := 

{!,... ,n}. For / G (^^^), let Ai(A) denote the maximal mi- 
nor of a k X n matrix A located in the column set /. The 
map A (A/(A)), where / ranges over (^^^), induces the 

Pliicker embedding Grkn ^ Rp(fe)~\ and the Ai{A) are 
called Pliicker coordinates. 

Definition 1. The totally non- negative Grassmannian {Grkn)>o 
(respectively, totally positive Grassmannian {Grkn)>o) is the 
subset of Grkn that can be represented by k x n matrices A 
with all A I (A) non-negative (respectively, positive). 

Postnikov [9 gave a decomposition of {Grkn)>o into 
positroid cells. For M C (t^^), the positroid cell S^J^ is the 
set of elements of {Grkn)>Q represented by all /c x n matrices 
A with the Ai(A) > for / G X and Aj(A) = 0, for J ^ M. 

Clearly {Grkn)>Q is a disjoint union of the positroid cells 
S^jOl" - in fact it is a CW complex p^. Note that {Grkn)>o 
is a positroid cell; it is the unique positroid cell in {Grkn)>Q 
of top dimension k(n — k). Postnikov showed that the cells of 
{Grkn)>Q are naturally labeled by (and in bijection with) the 
following combinatorial objects [9]: 

• Grassmann necklaces T of type (/c, n) 

• decorated permutations tt on n letters with k weak ex- 
cedances 

• equivalence classes of reduced plabic graphs of type (/c, n) 

• J-diagrams of type {k,n). 

For the purpose of studying solitons, we are interested 
only in the subset of positroid cells which are irreducible. 
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Definition 2. We say that a positroid cell Sm is irreducible if 
the reduced-row echelon matrix A of any point in the cell has 
the following properties: 

1. Each column of A contains at least one nonzero element. 

2. Each row of A contains at least one nonzero element in 
addition to the pivot. 

The irreducible positroid cells are indexed by: 

• irreducible Grassmann necklaces X of type (/c, n) 

• derangements tt on n letters with k excedances 

• equivalence classes of irreducible reduced plabic graphs of 
type (/c, n) 

• irreducible J-diagrams of type (k^n). 

We now review the definitions of these objects and some 
of the bijections among them. 

Definition 3. An irreducible Grassmann necklace of type (/c, n) 
is a sequence X — (/i, . . . , 7^) of subsets Ir of [n] of size k such 
that, for i G [n], /i+i = {li \ {i}) U {j} for some j / i. (Here 
indices i are taken modulo n.) 

Example 1. An example of a Grassmann necklace of type (4,9) 
is (1257, 2357, 3457, 4567, 5678, 6789, 1789, 1289, 1259). 
Definition 4. A derangement tt = (tti, . . . , vTn) is a permutation 
7T ^ Sn which has no fixed points. An excedance of tt is a pair 
(i^TVi) such that m > i. We call i the excedance position and 
TTi the excedance value. Similarly, a nonexcedance is a pair 
(i^TTi) such that tvi < i. 

Definition 5. A plabic graph is a planar undirected graph G 
drawn inside a disk with n boundary vertices 1, . . . , n placed in 
counterclockwise order around the boundary of the disk, such 
that each boundary vertex i is incident to a single edge\j Each 
internal vertex is colored black or white. 

Definition 6. Let Yx denote the Young diagram of the partition 
X. A d-diagram (or Le-diagram) L — (X,D)k,n of type (/c, n) 
is a Young diagram Yx contained in a kx (n — k) rectangle to- 
gether with a filling D : Yx ^ {0, +} which has the d-property: 
there is no which has a + above it in the same column and 
a -\- to its left in the same row. A d-diagram is irreducible if 
each row and each column contains at least one +. 
See Figure [1] for an example of an irreducible J-diagram. 
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k = 6, n= 16 

>^ = (10,9,9,8,5,2) 



Fig. 1. A Le-diagram (A, D)fc^^. 

Theorem 2. Theorem 17.2] Let S^j!^ be a positroid cell in 
{Grkn)>o- For 1 < r < n, let Ir be the index set of the mi- 
nor in M which is lexicographically minimal with respect to 
the order r<r + l < ••• <n< 1 < 2< ...r — 1. Then 
X[M) :— (/i, . . . , In) is a Grassmann necklace of type (/c, n). 
Lemma 1. /P, Lemma 16.2] Given an irreducible Grassmann 
necklace X, define a derangement tv = 7v{X) by requiring that: 
if li^i — [li \ {i}) U {j} for j 7^ i, then 7T{j) =2.0 Indices 
are taken modulo n. Then X 7r{X) is a bisection from irre- 
ducible Grassmann necklaces X — (/i, . . . , In) of type (/c, n) to 
derangements 7t{X) G Sn with k excedances. The excedances 
of Ti{X) are in positions Ii. 



Remark 1. // the positroid cell S^J^ is indexed by the Grass- 
mann necklace X, the derangement tv, and the d-diagram 
L, then we also refer to this cell as 5'i"",5'*"" and 6'!^^. 
The bijections above preserve the indexing of cells, that is, 

Qtnn Qtnn Qtnn 

— *^X(M) — *^7r(X(A4))' 

Soliton solutions to the KP equation 

Here we explain how to obtain a soliton solution to the KP 
equation from a point of {Grkn)>o- 

From the Grassmannian to the T -function. We start by fixing 
real parameters kj such that ki < k,2 • ■ ■ < Kn, which are 
generic, in the sense that the sums X^^^^ i^jm cire all distinct 
for 2 < d < /c. 

Let {Ej;j = 1, . . . , n} be a set of exponential functions in 
(x,2/,t) G defined by 

Ej{x, y, t) := exp [kjx + K^y + K^t) . 

If E^^-* denotes Ei/dx^ = i^iEi, then the Wronskian deter- 
minant with respect to x of , . . . , is defined by 

Wr(^i, ...,En) = det[{E^/~^^)i<ij<n] = H^^^'"^^) Ei ■ ■ ■ En 

i<j 

Let A be a full rank k x n matrix. We define a set of 
functions {/i , . . . , //e} by 

(/l,/2,...,A)^ = A-(^l,^2,...,^n)^, 

where (. . .)^ denotes the transpose of the vector (...)• The 
T -function of A is defined by 

rA{x,y,t) := Wr(/i,/2,...,/fc). [1] 

It's easy to verify that ta only depends on which point of 
{Grkn)>o the matrix A represents. 

Applying the Binet-Cauchy identity to the fact that fi = 
X^j=i ^ijEj for i = 1, . . . , /c, we get 

TA{x,y,t)= E Ai{A)Ei{x,y,t), [2] 

where Ei{x, y, t) with / = {kj^ , • • • , i^Jk } is defined by 

Ei := WriEj, , Ej„ . . . , E,,) = H -kj, ) Ej, ■ ■ ■ Ej, > 0. 

l<m 

Therefore if A G {Grkn)>o, then r a > for all {x,y,t) e M^. 

Thinking of ta as a function of A, we note from |l2j that 
the T function encodes the information of the Pliicker em- 
bedding. More specifically, if we identify each function Ei 
with / = {ji , . . . , j/c} with the wedge product Ej^ A • • • A Ej^^ , 

then the map r : Grkn ^ Rp(fe)~\ A ta has the Pliicker 
coordinates as coefficients. 

From the 7"-function to solutions of the KP equation. The 

KP equation 

dx \ dt dx dx^ J dy'^ 

was proposed by Kadomtsev and Petviashvili in 1970 [1], 
in order to study the stability of the one-soliton solution of 



""^The convention of 9 was to place the boundary vertices in clockwise order. 

^Actually Postnikov's convention was to set 7r(i) — j above, so the permutation we are associating 

is the inverse one to his. 



2 I www.pnas.org 



Footline Author 



the Korteweg-de Vries (KdV) equation under the influence of 
weak transverse perturbations. The KP equation also gives 
an excellent model to describe shallow water waves 

It is well known (see e.g. ^^at the r-function defined 

in [ij provides a soliton solution of the KP equation, 

UA{x,y,t) = 2—\nrA{x,y,t). [3] 
Note that if A G {Grkn)>o, then UA{x,y,t) is regular. 

From soliton solutions to soliton graphs 

One can visualize such a solution UA{x,y,t) in the x^-plane 
by drawing level sets of the solution for each time t. For each 
r G M, we denote the corresponding level set by 

Cr{t) :={{x,y)eR^ :uA{x,y,t)=r}. 

Figure [2] depicts both a three-dimensional image of a solution 
i^a(x, t), as well as multiple level sets Cr(0). Note that these 
levels sets are lines parallel to the line of the wave peak. 




Fig. 2. A line-soliton solution from A = (1, 1) G (Gri^2)>o- The left figure 
shows the 3-dimensional profile of UA{x,y,0)', the right one shows level sets of 
ua{x, 0). Ei represents the dominant exponential in the region. 

To study the behavior of ua{x, y, t) for A G 6*^^, we set 
fA{x,y,t) = max^{Aj{A)Ej{x,y,t)} 

= max I e'" (^^ (^)^^) +E ti -+4 ^+4 \ 
JeM\ J' 

where Kj := Yli<m('^jm ~ i^je) > 0- From |l2||, we see that 
generically, ta can be approximated by /a- 

Let fA{x,y,t) be the closely related function 

fA{x,y,t) = max <\n{Aj{A)Kj) + J^ii^Ji^ + i^hV + ^ii^) 1 • 
-^^^ L i^i ) 

[4] 

Clearly a given term dominates fA{x,y,t) if and only if its 
exponentiated version dominates fA{x,y,t). 
Definition 7. Given a solution ua{x, y, t) of the KP equation as 
in JSj; we define its contour plot CtoiuA) for each t = to to 
be the locus in where fA{x^y^t — to) is not linear. 
Remark 2. CtQ{uA) provides an approximation of the location 
of the wave crests. 

It follows from Definition [71 that CtQ{uA) is a one- 
dimensional piecewise linear subset of the x^-plane. 
Proposition 3. // each Ki is an integer, then CtoiuA) is a trop- 
ical curve in . 

Note that each region of the complement of Cto{uA) 
in R^ is a domain of linearity for /A(x,2/,to), and hence 
each region is naturally associated to a dominant exponential 
Aj{A)Ej{x,y,to) from the r-function |l2j. We call the line 
segments comprising CtoiuA) line-solitons^ Some of these 



line-solitons have finite length, while others are unbounded 
and extend in the y direction to ±oo. We call these unbounded 
line-solitons. Note that each line-soliton represents a balance 
between two dominant exponentials in the r-function. 

Lemma 2. pi Proposition 5] The dominant exponentials of the 
T -function in adjacent regions of the contour plot in the xy- 
plane are of the form E(i, 7712, . . . , ruk) and E(j, 7712, . . . , ruk)- 
According to Lemma O those two exponential terms have 
k — 1 common phases, so we call the soliton separating them 
a line-soliton of type [i,j]. Locally we have 

rA^Ai{A)Ei + Aj{A)Ej 

k 

= (Ai(A)KiE^ + Aj{A)KjEj)YlEm, 

with Ki = 11^=2 1^^ ~ I Tli<j l^^i ~ I' equation 
for this line-soliton is 

x+(Ki + Kj)y+{K^i+KiKj + K^j)t = — ^—In^^j^^. [5] 

Kj — Ki /\j[A)Kj 

Note that the ratio of the Pliicker coordinates labeling the 
regions separated by the line-soliton determines the location 
of the line-soliton. 

Remark 3. Gonsider a line-soliton given by [Sj. Gompute the 
angle ^[ij] between the line-soliton and the positive y-axis, 
measured in the counterclockwise direction, so that the nega- 
tive X-axis has an angle of ^ and the positive x-axis has an 
angle of Then tan ^[^^^j = Ki-\- kj . Therefore we refer to 
hvi + hvj as the slope of the [i,j] line-soliton (see Figure\^. 

We will be interested in the combinatorial structure of a 
contour plot, that is, the pattern of how line-solitons interact 
with each other. To this end, in Definition [9] we will associate 
a soliton graph to each contour plot. 

Generically we expect a point of a contour plot at which 
several line-solitons meet to have degree 3; we regard such 
a point as a trivalent vertex. Three line-solitons meeting at 
a trivalent vertex exhibit a resonant interaction (this corre- 
sponds to the balancing condition for a tropical curve). One 
may also have two line-solitons which cross over each other, 
forming an X-shape: we call this an X-crossing, but do not 
regard it as a vertex. In general, there exists a phase-shift 
at each X-crossing. However we ignore them in this paper as 
explained in the footnote 3. Vertices of degree greater than 4 
are also possible. 

Definition 8. A contour plot is called generic if all interactions 
of line-solitons are at trivalent vertices or are X -crossings. 

The following definition of soliton graph forgets the met- 
ric data of the contour plot, but preserves the data of how 
line-solitons interact and which exponentials are dominant. 

Definition 9. Let Cto{uA) be a generic contour plot with n un- 
bounded line-solitons. Golor a trivalent vertex black (respec- 
tively, white) if it has a unique edge extending downwards 
(respectively, upwards) from it. Label each region with the 
dominant exponential Ej and each edge (line-soliton) by the 
type [i,j] of that line-soliton. Preserve the topology of the 
metric graph, but forget the metric structure. Embed the re- 
sulting graph with bicolored vertices and X -crossings into a 
disk with n boundary vertices, replacing each unbounded line- 
soliton with an edge that ends at a boundary vertex. We call 
this labeled graph a soliton graph. 



"^In general, there exist phase-shifts which also appear as line segments (see 0). However the 
phase-shifts depend only on the k parameters, and we ignore them in this paper. 
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See Figure [3] for an example of a soliton graph. Although 
we have not labeled all regions or all edges, the remaining 
labels can be determined using Lemma O 




Fig. 3. Example of a soliton graph associated to S^^'^ with tt = 
(7, 4, 2, 8, 1, 3, 9, 6, 5). Each Eijj^i represents the dominant exponential in the 
r-function. This soliton graph was obtained from a contour plot by embedding it in 
a disk and coloring vertices appropriately. 

Permutations and soliton asymptotics 

Given a contour plot CtoiuA) where A belongs to an irre- 
ducible positroid cell and to is arbitrary, we show that the 
labels of the unbounded solitons allow us to determine which 
positroid cell A belongs to. Conversely, given A in the ir- 
reducible positroid cell S^'^ , we can predict the asymptotic 
behavior of the unbounded solitons in Cto{uA)- 
Theorem 4. Suppose A is an element of an irreducible positroid 
cell in {Grkn)>o- Consider the contour plot CtoiuA) for any 
time to. Then there are k unbounded line-solitons at y ^ 
which are labeled by pairs [cr^jr] with Cr < jr, and there are 
n — k unbounded line-solitons at y <^ which are labeled by 
pairs [ir^Qr] with ir < Qr- We obtain a derangement in Sn 
with k excedances by setting 7r(er) = jr and 7r(^r) = ir- More- 
over, A must be an element of the cell S^^^ . 

The first part of this theorem follows from work of 
Chakravarty and Kodama [5*, Prop. 2.6 and 2.9], ^ Theorem 
5] . Our contribution is that the derangement tt is precisely the 
derangement labeling the cell S^^^^ that A belongs to. This 
fact is the first step towards establishing that various other 
combinatorial objects in bijection with positroid cells (Grass- 
mann necklaces, plabic graphs) carry useful information about 
the corresponding soliton solutions. 

We now give a concrete algorithm for writing down the 
asymptotics of the soliton solutions of the KP equation. 
Theorem 5. Fix generic parameters ki < ■ ■ ■ < Kn- Let A be an 
element from an irreducible positroid cell <St^^ in {Grkn)>o- 
(So TT must have k excedances.) For any to, the asymptotic 




[2,4] [3,6] [4,8] 

Fig. 4. Unbounded line-solitons for tt = (6, 7, 1, 2, 8, 3, 9, 4, 5). Each Eij^i 
shows the dominant exponential in this region. 



behavior of the contour plot CtQ{uA) - i-e. its unbounded line- 
solitons, and the dominant exponentials in its unbounded re- 
gions - can be read off from tt as follows. 

• For y ^ 0, there is an unbounded line- soliton of type 
[z,7r(i)] for each excedance 7r(i) > i. From left to right, list 
these solitons in decreasing order of the quantity Ki + ^^^^(i) • 

• For y <^ 0, there is an unbounded line-solitons of type 
[7T{j),j] for each nonexcedance TT(j) < j. From left to right, 
list these solitons in increasing order of kj + k-kU) ■ 

• Label the unbounded region for x ^ with the exponential 
Ei^^,,,^ij^, where ii,. . . ,ik are the excedance positions of tt. 

• Use Lemma\^to label the remaining unbounded regions of 
the contour plot. 

Example 6. Consider the positroid cell corresponding to tt = 
(6, 7, 1, 2, 8, 3, 9, 4, 5) G Sg. The algorithm of Theorem \S\ gives 
rise to the picture in Figure 2} If one reads the dominant ex- 
ponentials in counterclockwise order, starting from the region 
at the left, then one recovers the Crassmann necklace X from 
Example{l\ Also note that tt{T) = tt. See Theorem^ 

Grassmann necklaces and soliton asymptotics 

One particularly nice class of positroid cells is the TP or totally 
positive Schubert cells. These are the positroid cells indexed 
by J-diagrams which are filled with all +'s, or equivalently, 
the positroid cells indexed by derangements tt such that tt"^ 
has at most one descent. When is a TP Schubert cell, we 
can make a link between the corresponding soliton solutions 
of the KP equation and Grassmann necklaces. 

Theorem 7. Let A be an element of a TP Schubert cell , 
and consider the contour plot Cto{uA) for an arbitrary time 
to . Let the index sets of the dominant exponentials of the un- 
bounded regions of Cto{uA) be denoted Ri,. . . ,Rn, where Ri 
labels the region at x <^ 0, and i?2, • • • , Rn label the regions in 
the counterclockwise direction from Ri. Then . . . , Rn) is 
a Crassmann necklace X and tt{X) — tt. 

Theorem [71 is illustrated in Example [6l 

Remark 4. Theorem^does not hold if we replace "TP Schubert 
ceir by "positroid cell." 



From soliton graphs to generalized plabic graphs 

In this section we associate a generalized plabic graph Pl{C) to 
each soliton graph C. We then show that from Pl{C) - whose 
only labels are on the boundary vertices - we can recover the 
labels of the line-solitons and dominant exponentials of C. 

Definition 10. A generalized plabic graph is a connected graph 
embedded in a disk with n boundary vertices labeled 1, . . . ,n 
placed in any order around the boundary of the disk, such that 
each boundary vertex i is incident to a single edge. Each inter- 
nal vertex must have degree at least two, and is colored black 
or white. Edges are allowed to form X -crossings (this is not 
considered to be a vertex). 

We now generalize the notion of trip from 9 , Section 13]. 

Definition 11. Civen a generalized plabic graph G, the trip T^ 

is the directed path which starts at the boundary vertex i, and 
follows the "rules of the road": it turns right at a black vertex, 
left at a white vertex, and goes straight through an X -crossing. 
Note that Ti will also end at a boundary vertex. The trip per- 
mutation TTG is the permutation such that ttq (i) = j whenever 
the trip starting at i ends at j . 

We use these trips to associate a canonical labeling of 
edges and regions to each generalized plabic graph. 
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Definition 12. Given a generalized plabic graph G with n bound- 
ary vertices, start at each boundary vertex i and label every 
edge along trip Ti with i. Such a trip divides the disk contain- 
ing G into two parts: the part to the left of Ti, and the part 
to the right. Place an i in every region which is to the left of 
Ti. After repeating this procedure for each boundary vertex, 
each edge will be labeled by up to two numbers (between 1 and 
n ), and each region will be labeled by a collection of numbers. 
Two regions separated by an edge labeled ij will have region 
labels S and (S \ {i}) U {j}. When an edge is assigned two 
numbers i < j, we write [i,j] on that edge, or {i,j} or {j,i} 
if we do not wish to specify the order of i and j . 
Definition 13. Fix an irreducible cell of {Grkn)>o- To 

each soliton graph G coming from a point of that cell we as- 
sociate a generalized plabic graph Pl{G) by: 

• labeling the boundary vertex incident to the edge {i,7Vi} by 

7Vi = 7v{i), 

• forgetting the labels of all edges and regions. 

See Figure [5] for the generalized plabic graph Pl{C) corre- 
sponding to the soliton graph G from Figure [3] 

Theorem 8. Fix an irreducible cell S^^^ of {Grkn)>o, and con- 
sider a soliton graph G coming from a point of that cell. Then 
the trip permutation associated to the plabic graph Pl{G) is 
TV, and by labeling edges and regions of Pl{G) according to 
Definition \12\ we will recover the original labels in G . 




Fig. 5. Example of a generalized plabic graph G{C). 

We invite the reader to apply Definition [12] to Figure O and 
then compare the result to Figure [S] 

Remark 5. By Theorem we can identify each soliton graph 
G with its generalized plabic graph Pl{G). 

Soliton graphs for positroid cells when t <^ 

In this section we give an algorithm for producing a general- 
ized plabic graph G- (L) from the J-diagram L of a positroid 
cell S^i^^ . It turns out that this generalized plabic graph gives 
rise to the soliton graph for a generic point of the cell Sj^, 
at time t ^ sufficiently small. 

Algorithm 9. Given a d-diagram L, construct G-{L) as fol- 
lows: 

1. Start with a J-diagram L contained in a k x {n — k) rect- 
angle. Label its southeast border by the numbers 1 to n, 
starting from the northeast corner. Replace ^s and + 's by 
"crosses^^ and "elbows^\ From each label i on the southeast 
border, follow the associated "pipe^^ northwest, and label its 
destination by i as well. 



2. Add an edge, and one white and one black vertex to each 
elbow, as shown in the upper right of Figure Forget the 
labels of the southeast border. If there is an endpoint of a 
pipe on the east or south border whose pipe starts by going 
straight, then erase the straight portion preceding the first 
elbow. 

3. Forget any degree 2 vertices, and forget any edges of the 
graph which end at the southeast border of the diagram. 
Denote the resulting graph G-{L). 

4. After embedding the graph in a disk with n boundary ver- 
tices, we obtain a generalized plabic graph, which we also 
denote G-{L). If desired, stretch and rotate G-{L) so that 
the boundary vertices at the west side of the diagram are at 
the north instead. 
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Fig. 6. AlgorithmlUfor S^j^"^ with L the Le-di agram from the upper left. 



Figure [6] illustrates the steps of Algorithm [9] Note that 
this produces the graph from Figure [5] 

Theorem 10. Let L be a d-diagram and tv = 7v{L). Then 
G-{L) has trip permutation tt. Label its edges and regions 
according to the rules of the road. When S^^^ is a TP Schu- 
bert cell, then G-{L) coincides with the soliton graph Gt{uA), 
provided that A G S^^^ is generic and t <^ sufficiently small. 
When S^L^ is an arbitrary positroid cell, we can realize G-{L) 
as ^^most" of a soliton graph Gt{uA) for A G S^^^ and t <^ 0. 
Moreover, we can construct Gt{uA) from G-{L) by extending 
the unbounded edges of G-{L) and introducing X -crossings as 
necessary so as to satisfy the conditions of Theorem\^ 

Reduced plabic graphs and cluster algebras 

The most important plabic graphs are those which are re- 
duced [9 , Section 12]. Although it is not easy to characterize 
reduced plabic graphs (they are defined to be plabic graphs 
whose move- equivalence class contains no graph to which one 
can apply a reduction)^ they are important because of their 
application to cluster algebras and parameterizations of cells. 
Theorem 11. Let A be a point of a TP Schubert cell, let to be 
an arbitrary time, and suppose that the contour plot Cto{uA) 
is generic and has no X-crossings. Then the soliton graph 
associated to Cto{uA) is a reduced plabic graph. 

Cluster algebras are a class of commutative rings with a 
remarkable combinatorial structure, which were defined by 
Fomin and Zelevinsky [13 . Scott [M] proved that Grassman- 
nians have a cluster algebra structure. 

Theorem 12. [T^ The coordinate ring of the (affine cone over 
the) Grassmannian has the structure of a cluster algebra. 
Moreover, the set of labels of the regions of any reduced plabic 
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graph for the TP (totally positive) Grassmannian comprises a 
cluster for this cluster algebra. 

Remark 6. Scott's strategy in Jl^J was to show that certain la- 
belings of alternating strand diagrams for the TP Grassman- 
nian gave rise to clusters. However, alternating strand di- 
agrams are in bijection with reduced plabic graphs 19], and 
under this bijection, Scott 's labelings of alternating strand di- 
agrams correspond to the labelings of regions of plabic graphs 
induced by the various trips in the plabic graph. 

Corollary 13. The Pliicker coordinates labeling regions of a 
generic soliton graph with no X-crossings for {Grk,n)>o form 
a cluster for the cluster algebra o/G'rfc,n. 

Conjecturally, every positroid cell <S*^^ of the totally non- 
negative Grassmannian also carries a cluster algebra struc- 
ture, and the Pliicker coordinates labeling the regions of any 
reduced plabic graph for S^^^ should be a cluster for that 
cluster algebra. In particular, the TP Schubert cells should 
carry cluster algebra structures. Therefore we conjecture that 
Corollary [13] holds with "TP Schubert cell" replacing "TP 
Grassmannian." Finally, there should be a suitable general- 
ization of Corollarv 113! for arbitrary positroid cells. 



The inverse problem 

The inverse problem for soliton solutions of the KP equation is 
the following: given a time t together with the contour plot of 
a soliton solution, can one reconstruct the point of {Grk,n)>o 
which gave rise to the solution? 

Theorem 14. Fix ki < ■ ■ ■ < as usual. Gonsider a generic 
contour plot of a soliton solution coming from a point A of 
a positroid cell , for t ^ 0. Then from the contour plot 
together with t we can uniquely reconstruct the point A. 

The strategy of the proof is as follows. From the contour 
plot together with t, we can reconstruct the value of each of 
the dominant exponentials (Pliicker coordinates) labeling re- 
gions of the graph. We have shown how to use the J-diagram 
to construct the soliton graph for a positroid cell when t <C 
is sufficiently small, which allows us to identify what is the set 
of Pliicker coordinates which label regions of the graph. We 
then show that this collection of Pliicker coordinates contains 
a subset of Pliicker coordinates, which Talaska 15j showed 
were sufficient for reconstructing the original point of S^^'^. 

Using Theorem [Ml Corollary [131 and the cluster algebra 
structure for Grassmannians, we can solve the inverse problem 
for the TP Grassmannian for any time t. 
Theorem 15. Gonsider a generic contour plot of a soliton so- 
lution coming from a point A of the TP Grassmannian, at 
an arbitrary time t. If the contour plot has no X-crossings, 
then from the contour plot together with t we can uniquely 
reconstruct the point A. 
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Triangulations of a polygon and soliton graphs 

We now explain how to use triangulations of an n-gon to pro- 
duce all soliton graphs for the TP Grassmannian {Gr2,n)>o- 

Algorithm 16. Let T be a triangulation of an n-gon P, whose 
n vertices are labeled by the numbers 1, 2, ... ,71; in counter- 
clockwise order. Therefore each edge of P and each diagonal 
of T is specified by a pair of distinct integers between 1 and 
n. The following procedure yields a labeled graph ^(T). 

1. Put a black vertex in the interior of each triangle in T . 

2. Put a white vertex at each of the n vertices of P which 
is incident to a diagonal of T; put a black vertex at the 
remaining vertices of P. 

3. Gonnect each vertex which is inside a triangle of T to the 
three vertices of that triangle. 

4. Erase the edges of T , and contract every pair of adjacent 
vertices which have the same color. This produces a new 
graph G with n boundary vertices, in bijection with the ver- 
tices of the original n-gon P. 

5. Add one unbounded ray to each of the boundary vertices 
of G, so as to produce a new (planar) graph ^(T). Note 
that ^(T) divides the plane into regions; the bounded re- 
gions correspond to the diagonals of T , and the unbounded 
regions correspond to the edges of P. 

6. Resolve any non-trivalent vertices into trivalent vertices. 

Theorem 17. The graphs ^(T) constructed above are soliton 
graphs for {Gr2,n)>o- Gonversely, any generic soliton graph 
with no X-crossings for {Gr2,n)>o comes from AlaorithmHR 




Fig. 7. Algorithmic starting from a triangulation of a hexagon. 



Flipping a diagonal in a triangulation corresponds to a 
mutation in the cluster algebra. In our setting, each muta- 
tion may be considered as an evolution along a flow of the KP 
hierarchy defined by the symmetries of the KP equation. 
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